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We scrutinize the diagrammatic perturbation theory of noninteracting electrons in a random 
potential with the aim to accomplish a consistent comprehensive theory of quantum diffusion. Ward 
identity between the one-electron self-energy and the two-particle irreducible vertex is generally 
not guaranteed in the perturbation theory with only elastic scatterings. We show how the Ward 
identity can be established in practical approximations and how the functions from the perturbation 
expansion should be used to obtain a fully consistent conserving theory. We derive the low-energy 
asymptotics of the conserving full two-particle vertex from which we find an exact representation of 
the diffusion pole and of the static diffusion constant in terms of Green functions of the perturbation 
expansion. We illustrate the construction on the leading vertex corrections to the mean-field diffusion 
due to maximally-crossed diagrams responsible for weak localization. 
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I. INTRODUCTION 

Diffusion of particles drifted in inhomogeneous materi¬ 
nis by external forces manifests itself on large distances 
and hence it is a macroscopic phenomenon. That is why 
it would be treated mostly classically by means of the 
diffusion equation or variants of the Boltzmann trans- 
port equation.lil Observed deviations from the behavior 
predicted by the Boltzmann equation are of quantum 
origin. To explain them we have to develop a microscopic 
quantum theory of diffusion where elementary objects 
are waves. Their diffusion is caused by wave scatterings 
on irregularly distributed impurities in otherwise regular 
environment. Possible interference of waves makes the 
qualitative difference between the classical and quantum 
diffusion. A theory of multiple scatterin gs o f waves on a 
random potential was introduced by Lax.l212l The theoreti- 
cal framework for the microscopic description of quantum 
transport and diffusion was set by the linear-response 
theory of Kubo.® The interest in microscopic diffusion of 
quantum particles in disordered Systems was boosted by 
P. W. Anderson who discovered that charge diffusion due 
to multiple electron scatterings on randomly distributed 
impurities can vanish if the strength of the random po¬ 
tential is sufhciently large.^^ The full understanding of the 
Anderson localization transition, that is, the tr ans ition 
from diffusive to insulating phase, is stili elusive.l^ 

There are two major streams in the study of suppression 
and eventual v anis hing of diffusion: finite-cluster numeri- 
cal simulationi^ and (semi-)analytic approaches in the 
thermodynamic limit.l^Unfortunately, agreement between 
the results and conclusions of the numerical simulations 
and the analytic theories are not yet fully satisfactory. 
Each approach has its advantages and drawbacks. The 
numerical simulations give rather accurate values of the 
boundary of the localized phase and of the critical ex- 
ponents, their results need, however, to be appropriately 
scaled to get rid of finite-size effects. The theoretical 
approaches in the thermodynamic limit directly access 


thermodynamically relevant quantities, but suffer from 
the fact that only homogeneous quantities are available 
and some of the relevant questions cannot be addressed. 

So far, most of the theoretical approaches to quantum 
transport work well either in the metallic regime, far 
from the localization transition, or just in the critical 
region of the transition. Perturbation, diagrammatic ex- 
pansions in powers of the random potential formulated in 
the language of the averaged Green functions have been 
used to describe the metallic phase.^ Transformations to 
effective solvable models then address long-range ffuctua- 
tions and the critical behavior of correlation and response 
functions at the Anderson localization transition in the 
strong-disorder limit.lS^ Neither of these approaches can 
smoothly be extended to the opposite limit. The diagram¬ 
matic approach works well for one-electron functions and 
for equilibrium thermodynamic properties but even its 
self-consistent extension fails to incorporate adequately 
vertex corrections to the semi classi cal Drude expression 
for the electrical conductivity.li2*i^ The latter construc¬ 
tion does not cover adequately the non-universal and 
non-critical properties of the disordered Systems away 
from the localization transition.^ It is the diagrammatic 
approach that is fully microscopic and has the potential to 
offer a fully consistent first-principles theory of quantum 
diffusion. 

There are presently two diagrammatic approaches to 
the microscopic description of quantum diffusion aim- 
ing at covering both metallic and insulating phases. It 
is the self-consistent theory of Anderson localization of 
Vollhardt and Wolfl^^^ and the parquet approach of 
the present authors.^ ^^h^ l xhe common quality of both 
approaches is a two-particle self-consistency used to de¬ 
scribe the Anderson localization transition.^Sl They differ, 
however, in the way this self-consistency is used. The for¬ 
mer approach uses the framework delimited by the Ward 
identity and introduces a self-consistent equation for the 
diffusion constant governing the long-range fluctuations 
of the singular part of the electron-hole correlation fune- 
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tion. Macroscopic conservation laws are obeyed in this 
approach. A systematic and controllable way of improv- 
ing the self-consistent equation for the diffusion constant 
is, however, missing. The latter approach allows one to 
develop controllable approximations for the two-particle 
vertices including two-particle self-consistency, but it is 
unable to reconcile the calculated vertex with the Ward 
identity and macroscopic conservation laws. 

The aim of this paper is to present a systematic and 
analytically controllable diagrammatic description of dif- 
fusive transport in models of non-interacting electrons 
scattered on randomly distributed impurities. The non- 
local two-particle scattering processes calculated in the 
full two-particle state space lead to vertices and Green 
functions that generically do not comply with the Ward 
identity and hence do not constitute a conserving theory.l^ 
To arrive at a conserving approximation, the two-particle 
irreducible vertex from the diagrammatic construction is 
only an auxiliary function from which the physical one 
is appropriately constructed. We show how to make di¬ 
agrammatic theories conserving and how to restore the 
Ward identity in approximate treatments. The physical 
irreducible vertex obeying the Ward identity is then used 
in a Bethe-Salpeter equation to determine the physical 
two-particle functions from which all measurable quan- 
tities are determined. Thereby, the construction of the 
two-particle functions becomes fully consistent and the 
relevant macroscopic thermodynamic relations in the dif- 
fusive regime hold.^ We further derive an exact form 
of the diffusion pole and of the diffusion constant and 
apply them to derive leading non-local corrections to the 
mean-field diffusion constant. 

The paper is organized as follows. We introduce the 
model, notation and basic relations in Sec. II. The way the 
dynamical Ward identity is established in perturbation 
theory is presented in Sec. III. We derive in Sec. IV a low- 
energy singularity in the full two-particle vertex and with 
the aid of the dynamical Ward identity also the diffusion 
pole in the electron-hole correlation function. The exact 
form of the static diffusion constant and its approximate 
form including maximally crossed diagrams are presented 
in Sec. V. Section VI brings concluding discussion. 


II. DEFINITIONS, NOTATION AND BASIC 
RELATIONS 

We will consider a noninteracting lattice electron gas 
scattered on random impurities described by the Anderson 
tight-binding Hamiltonian 

H = 'y ^ tijc^cj + 'y ^ i^cjci, (I) 

<ij> i 

where Vi is a local, site-independent random potential. 
We assume validity of the ergodic theorem leading to 
the existence of the thermodynamic limit, in which all 
quantities are homogeneous, configurationally averaged. 
We use the Standard diagrammatic representation of the 


scattering processes.l^ The state space is a Hilbert space 
spanned over the Bloch waves '0k(A) with energy E and 
wave vectors k from the first Brillouin zone. 

The fundamental function carrying all Information 
about the linear response of the disordered electron system 
is the two-particle Green function 
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( 2 ) 


where the brackets ()av indicate the configurational aver- 
aging. The indices i,j,... correspond to lattice sites with 
positions Ri, R^-.... Since the averaging recovers homo- 
geneity we use the Fourier transform to the momentum 
space to label the averaged Green functions. 

Only real energies are relevant for physical measurable 
quantities. Nevertheless, we have to distinguish the way 
the real energy is reached from the complex plane to 
obtain unambiguous results. Energies with an infinitesimal 
imaginary part 77 will be tacitly assumed in all expressions, 
and superscript R/A will correspond to positive/negative 
rj. We will use the following notation for the Fourier 
transform of the two-particle Green function {h = 1) 


ijkl 

X gi(k'+q/2).Rjg-i(k'-q/2)-Rfcgi(k-q/2)-Ri 

X Gi%i(E + w/2 + iO+,E- a;/2 - i0+). (3) 

The independent incoming and outgoing variables of the 
electron are k+ = k-|-q/2 and k^ = k'-|-q/2, and q is the 
difference between the incoming (outgoing) momentum 
of the electron and the hole. 

The Principal function of interest in studying quantum 
diffusion is the electron-hole correlation function. It dis- 
plays a low-energy singularity, a diffusion pole controlling 
the long-range electron diffusiori^ 

$(A;a;,q) = ^^GS^i(A;a;,q) 

kk' 

_, __ (4) 

tj-s-o.ij-s-o —iuj + D{Lo)q^ ’ 

where we denoted np the density of States at the Fermi 
energy E and D{ui) the dynamical diffusion constant. The 
electron-hole correlation function contains, however, only 
a reduced Information and is not an integral part of the 
perturbation theory. It is the two-particle Green function 
G^^,{E; oj, q) that comes out of the two-particle perturba¬ 
tion (diagrammatic) theory. It is the primary objective of 
the perturbation theory of quantum diffusion to identify 
the scattering events contributing to G^^,{E;u},c[) rele¬ 
vant for recovering the singular low-energy asymptotics 
of the electron-hole correlation function $(i?;a;,q). 

To determine the two-particle Green function one has 
to sum contributions from multiple scatterings to the 
one- and two-particle functions. To do so effectively one 
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introduces the concept of irreducibility that is best de- 
fined diagrammatically. The one-particle irreducibility is 
uniquely defined and means that cutting a single one- 
electron propagator does not disconnect the diagram. 
The contribution from all one-particle irreducible dia- 
grams is contained in the self-energy that renor- 

malizes the one-electron propagator via the Dyson equa- 
tion. The two-particle irreducibility is, however, ambigu- 
ous and de mands d istinction between local and nonlocal 
scatteringsfi^ESllIl Here we do not go into these subtle 
details of the two-particle perturbation theory. We intro- 
duce only the vertex functions that are necessary for the 
determination of the electron-hole correlation function. 
We first single out the uncorrelated propagation from the 
two-particle Green function G^^(E;u!,q) and introduce 
a two-particle vertex w, q), 

(E; o;, q) = (E+)Gt {E-) [fV4k' 

+ [E- cc, q)G«^ {E+)Gti._ {E.)] , (5) 

where we shortened the notation with the aid of defi- 
nitions k± = k ± q/2 and E± = E E uj/2. The vertex 
r^;^jj;^,(i?;w,q) contains only the correlated propagation. 
It can further be split into reducible and irreducible contri- 
butions from simultaneous scatterings of pairs of particles. 
The electron-hole irreducibility is introduced in the dia- 
grammatic representation so that a two-particle diagram 
cannot be disconnected by cutting just a single pair of 
electron-hole (antiparallel) lines. The full vertex is a sum 
of the irreducible and reducible contributions that can be 
represented by a Bethe-Salpeter equation 

(i?;a;,q) 

= {E- 0^, q) + ^ E 

xG^,,{E+)Gt,,{E.)T^f^,^[E-uj,q) . ( 6 ) 

Function L^^{E, uj; q) is the electron-hole irreducible 
vertex containing scatterings that are connected by more 
than One pair of electron-hole lines. This two-particle 
function is not independent of the one-particle self-energy. 
Vollhardt and Wolfl^i^ proved that if a Ward identity 

AE^^iE;u;,q) 

= ^ E q) (7) 

holds then macroscopic conservation laws are obeyed. Here 
we introduced AG^^(E-u;,q) = G^ (G+l - Gj^ (E_) 
and AE^^{E;u;,q) = E^JE+) - E^_iE_). 

Perturbation theory of quantum diffusion has to address 
simultaneously both one- and two-particle irreducible func¬ 
tions, and We demonstrated in our earlier 

publications that it is generally impossible to guarantee 
validity of the dynamical Ward identity in the pertur¬ 
bation expansion if we go beyond the local mean-field 


approximation of the irreducible functi ons, sinc e it is in 

conflict with causality of the self-energy! ^^b8 | 22 ] means 

that Eq. 0 is generally not obeyed in its full extent. A 
digression from the Ward identity is a universal feature of 
any expansion around the local mean-field solution. The 
problem is that the Ward identity dynamically restricts 
the correlated movement of two particles in an uncontrol- 
lable way.^ Ward identity in noninteracting system is not 
a consequence of a microscopic symmetry transformation, 
since the one- and two-particle functions are not dynami¬ 
cally coupled. There is hence no microscopic law forcing 
the Ward identity in random noninteracting Systems be 
strictly obeyed in virtual processes represented by individ- 
ual scattering events. To make the perturbation expansion 
accomplishable we then relax the full dynamical Ward 
identity when summing diagrammatic contributions to the 
vertex functions. We denote A^"^ the irreducible vertex 
that is obtained from the perturbation expansion and F^^ 
the full two-particle vertex resulting from the correspond- 
ing Bethe-Salpeter equation. We must strictly distinguish 
the vertex functions used in the two-particle perturbation 
theory from the thermodynamically consistent vertices 
that do not contradict macroscopic conservation laws. 


III. ESTABLISHING WARD IDENTITY IN THE 
PERTURBATION THEORY 

The Ward identity is a microscopic condition for making 
the theory conserving. Quantum theory contains virtual 
processes represented by Feynman diagrams for which 
macroscopic conservation laws need not hold unless forced 
by microscopic (local) symmetry transformations. There 
are no such symmetries in noninteracting Systems. It is 
only important and necessary that the conservation laws 
are restored in measurable quantities. To achieve this we 
reconcile the one-electron Green functions with the two- 
particle vertex via the Ward identity in the best possible 
way. The full dynamical Ward identity can neither be 
used to determine the one-particle self-energy from the 
two-particle irreducible vertex nor vice versa, since the 
vertex contains more information than the self-energy. 
The Ward identity generally serves only as a consistency 
check and a guarantee that the macroscopic conservation 
laws are obeyed. Ward identity Q for w = 0 and q = 0 
can nevertheless be used to determine the imaginary part 
of the self-energy from the electron-hole irreducible vertex 

QE^{E) = 1^AS(5(E;0,0)5G«(F;), ( 8 a) 

k' 

since both sides of this identity contain the same number 
of degrees of freedom and the equation can consistently 
be resolved for each energy E and momentum k. The 
corresponding real part is then found from the Kramers- 






Kronig relation 


IV. SINGULAR BEHAVIOR OF 
TWO-PARTICLE FUNCTIONS 


4 


/ OO 

-OO 


du! 92 ^( 0 ;) 
TT LU — E 


(8b) 


A. Low-energy asymptotics of the full two-particle 
vertex 


that ensures analyticity and causality of the self-energy 
in the plane of complex energies beyond the real axis. 

Using the one-electron Green functions with the self- 
energy from Eqs. ([^ in the equations for the two-particle 
vertices is not enough to guarantee validity of the full 
Ward identity. The Ward identity for nonzero frequen- 
cies and transfer momenta must be obeyed to guarantee 
macroscopic conservation laws. The non-local two-particle 
irreducible vertex Ak^/ (E; w, q) does not obey the full 
Ward identity and hence it has to be appropriately modi- 
fied on a subspace on which its action is determined by 
the Ward identity. We introduce a new function measur- 
ing the deviation of the given vertex from the Ward 
identity 


RuiE; q) = ^ X! q) AGk' (E; w, q) 

k' 

- AEk(i^;w,q) . (9) 

This function vanishes in the metallic phase for w = 0 
and <7 = 0 due to the definition of the self-energy, Eq. (|^ , 
used in the one-electron Green functions determining the 
vertex With the aid of function i?k we construet a 
new electron-hole irreducible vertex 


tRA 

^k+k; 


(£;;w,q) = A 


RA 

k+k 


+ 


(E;w,q) 


1 

(AG(E;..,q)2) 


X 


AGk(E;a;,q)i?k'(^^;w,q) -f i?k(^^; w, q)AGk'(E; w, q) 


AGk(E; w, q) AGk' (G; w, q) 
(AG(£;;cc,q)2) 


(i?(E;a;,q)AG(E;a;,q)) 

( 10 ) 


Validity of the dynamical Ward identity, Eq. Q’ is 
essential for the existence of a low-energy singularity in the 
two-particle vertex. Singularity of vertex (i?;a;,q) 

has, however, a more complicated structure than the 
diffusion pole in the electron-hole correlation function in 
Eq. Q, since the fermionic momenta k, k' remain relevant 
and we cannot disregard them. 

To derive the exact singular asymptotics of the full two- 
particle vertex in the metallic phase we assume analyticity 
of the irreducible vertex (E; w, q) in the low-energy 

limit w —>■ 0 and <7 —>■ 0 so that an expansion in powers of 
uj and q exists. As a first step we rewrite the product of 
two Green functions as a fraction containing differences 
of one-particle functions, 

g«X'B+)g1(e_) 

^_ AGk(E;a;,q) _ 

ASk(E; uj,q)-uj + Aek(q) ' 

We next use a small-momentum expansion 


/k±q/ 2 (z) = fk{z) ± i(7(q • Vk).fk(z) 

+ ^9^(q-Vk)(q-Vk)/k( 2 :) + 0((7^) (13) 

to get 


GtiE+)GtiE.) = 


AGk(E;a;,q) 


ASk(E;w,q) _ 
w - Aek(q) Aek(q)^ 


ASk(E;0,q) ASk(E; 0,0)2 


where 


(14a) 


that now obeys Eq. Q. We abbreviated (AG(E; w, q)^) = 
A^”^Ek AGk(E;a;,q)2 and {R{E;uj,q)AG{E;uj,q)) = 
N~^J2kE^(E;u},q)AGkiE]Uj,q). Function is the 
desired physical irreducible vertex to be used in Eq. ([^ to 
determine the physical value from which all relevant 
macroscopic quantities will be calculated. 

Since the Ward identity cannot be fully obeyed in the 
diagrammatic perturbation theory, vertex functions T and 
T differ and are equal only when the difference function 
Rk{E;ui,q) vanishes. It happens for w = 0 and q = 0, 
that is. 


r«k^(ii;;0,o) = r«k^(ii;;0,o), (ii) 

assuming the self-energy is determined from the two- 
particle vertex via Eqs. Q. It means that the vertex 
function A^^,{E-,oj,q) is directly related to the measur- 
able macroscopic quantities only for w = 0 , q = 0. 


Aek(q) = g(q-Vk) -l-O(g^) 
AEk(E; UJ, q) = 2iQY,^{E) + 


(14b) 


+ 9(q-Vk)3?E« + 0(g'), (14c) 


and 


AGk(E; UJ, q) = 2iQG^{E) + uj^G^ 

+ (7[(3?G«)2 - (5G^)2]q • (vk + Vk3?E^) 

-2 3?G^5G^(q. Vk)3E« + 0(g'). (14d) 

The dot indicates the derivative with respect to uj. We 
further introduce two new vertex functions 

(B; q) = (E; 0 ., q) . ( 15 ) 
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for which we formulate the Bethe-Salpeter equation in 
the low-energy limit. The leading low-energy asymptotics 


nRA 

y\ 


k+k'^ 




_ rRA 

— -^kk" 


-E 

N ^ 


rRA 

•^kk" 




Aek»(q) 

AEk"(i?;0,q) 


A£k»(q)^ ' 
AE2„ _ 


+ gq • (q) 


A£k''(q) \ 
AEu,, ) 


+ 'V(q-v,)^A 




5^4, {E-u:,q). 


(17) 


The functions without the values of their variables are taken at the Fermi energy E for w = lO"*" and q = 0. Moreover, 
momentum q at the irreducible vertices is set zero after the q derivatives were performed. After regrouping the terms 
we obtain 




k" 


■ W 


rRA 

•^kk" 


1 


AEw» 


rRA 

•^kk" 


q ■ Vk" 

AEk^' 


q • V, 


rRA 

‘-'k+K 


(q) 


q ■ Vk" / q ■ Vk" - q ■ VgSftsg, 

AEk" V AEk" 


^(q-v,)^ 


rRA 

-^k+k" 



G 


RA 

k"k 


+ 


(F;;u;,q) 


(18) 


We now split the two-particle vertex into two parts, one with even and one with odd symmetry with respect to 
momentum inversion q -A —q. We gather the odd contributione to form even-symmetry terms that stay in the 
denominator of both functions. Since we are interested only in the leading non-vanishing order we sum the contributione 
to the denominator of the full vertex to second order. The corresponding expansion for the even part is 


^^kEv(^;^>q)-^kk^' + AfE]^: 


N ■ 


-RA 

-kk" 


rRA ^ I rRA 

AEk-- 


J ^2 2-^ '^k+ki+lqj q q 


kik2 


AEu 


1 - . A^lk'' (q) (V, • q - V 


I kik2 ' 


-g2£«^4,,(q) 


^ , q • Vk" - q • 


VSfSjJ,, 


q • Vk" l,t ...2 


AE, 


(£;;w,q), (19) 


while that for the odd one reads 

^k+k' q) = ^ E (q) • q - 


q ■ Vk" 

AEk" 


rRA 

»k"k' 




-RA 

"kk" 


• (jJ 


+£• 


-RA 

kk" 


^ E ^kE+(q) -q- [l-^"^^(w) 


kik2 


-q^L^%,,(q) 


V, • q- 


AEk, 

q • Vk" - q • V3?E«, 

AEk" 


k 

rRA 

Jkik2 


rRA ^ 


(q) (V, 


q 


q • Vk' 


q • Vk" 1 _ 

AEk" 2^ 


AEk" 

2'Ek"k'(^^;w,q)- (20) 


■+ 


r 


The inverse of operator 1 — should be taken as the 
limit lim(^_j.o[ 1 — £^'^(w)] in order to regularize its 
low-energy pole. The same holds also for vertex Q^^{uj). 
The term containing the inversion [l — >C^"^(a;)] is 
dominant in the small frequency and small momentum 
expansion of the irreducible vertex £k+k' iE;uj,q) when 
applied in low spatial dimensions. 

The asymptotic expressions will now be used to deter- 
mine the low-energy behavior of the electron-hole correla- 


tion function and to find exact representations of the the 
diffusion pole and the diffusion constant. 


B. Low-frequency limit of the homogeneous 
two-particle vertex 

We first set q = 0 and investigate the low-frequency 
limit of the full vertex rkk'(w)- The homogeneous two- 
particle vertex is then determined from a Bethe-Salpeter 
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equation 

rkk'(w) = ivkk'(^) 

1 Lkk" (w)AGk" (w) 

AEk"(a;)—w 


rk"k'(w), (21) 


where we skipped the superscript RA at the two-particle 
vertices. It is straightforward to derive the following low- 
frequency asymptotics 


Ekk^(a^) — Ekk' “t“ tuJ^Ekk^ ; 
QLkk' = QAkk' 


2(5G2) L 


+i?kSGk' - (i?9G) 


9Gki?k' 

9GkSGk' 


\ / (3G2) 


(22a) 


(22b) 


k' 

+ ^E^kk'3?Gk--3?Sk 

3?Gk' 


AGk(a;) = 2i5Gk 

AEk(a;) = 2i3Ek 
AGk"(w) 3Gk 


AEk"(a;) — w 3Ek 


1 + w 

1 + w 

1 + 


5Gk 

SREk ' 

9Sk 

LO 


2i3Ek 




(22c) 

(22d) 

(22e) 

(22f) 


Using these results we expand Eq. (211 up to the leading 
small-frequency term 

rkk'(w) = Ekk' + 1^7 Akk» rk"k' (a^) 


3Ek^, 


N 


E 


• L\c\c" ( ^Ek'' ^ ^ • 


X |f^rk^'k'(a;). (23) 

AEk” 


We skipped the variables if their values were zero and the 
functions were evaluated at the Fermi energy for w = i0“'" 
and q = 0. 

The operator (matrix) form of the Bethe-Salpeter equa¬ 
tion for vertex Qiuj) reads 


Q{io) — G + C,Q -\- iuj 

with 

QGkk' 


3£-G 


1 


29E 


g{u:) (24) 


9Lkk' - 


Gkk' ( 5iGk' JiSi 


9Gk' 9E. 


SGk> 

3£k' 


(25) 


We demonstrate the existence of a pole in a specific ma¬ 
trix element of vertex Q(oj) in the limit w —>■ 0. We 
multiply Eq. ( |M| by a left normalized vector (9G| = 
{Nyj (S5G2))“^ S5Gk(k| and sum over the right mo¬ 
mentum in the two-particle vertices. If we introduce 
|^(w)) = A“2 X)k.k'|J^)^kk'(w) we obtain 


(3G|e(cc)) = 


(^G) 


(3G|e(a;)) 


+ iuj { 3G 


3£-G 


1 


23E 


0 (cc)) . (26) 


We assume that (5G| is the only maximal left eigenvector 
of operator C. We choose an appropriate orthonormal 
basis with the maximal eigenvector in the one-particle 
Hilbert space, that is, we have a decomposition of the 
unity operator 1 = |3G)(3G|+^„ |^„) {4>n\- We evaluate 
projections of vertex Giuj) to the basis vectors \(j)n) and 
obtain 


(<(>n|l/(w)) = {4>n\J^) + 
3G-G ^ 


£ 


Gico) 


25E 


GH) = 


£ 


5G 


X (5G|g(a;)) + ^ (</.„ |£| (</>m|0(w)) . (27) 

m 

We could neglect the first term on the right-hand side 
and, since we assumed that the operator of the frequency 

derivative 9£ is bounded, also the frequency-dependent 
term in the second line. They are nonsingular. The ab¬ 
solute term may be divergent via a projection onto the 
eigenvector |3G). Its divergent part reads: 


(<^„| 0 (a;))=^E„ (T-Pg) 


X (cj). 


9G) (3G|e(a;)) (28) 


£ 


where we denoted projectors P = |(/)„) ((^„| = l — Q = 

1 — |5G)(3G| and used an eigenvalue equation QC = Q 
together with the orthogonality relations ((/>„|3G) = 0. 

After summing over the intermediate States and using 
Eq. (27) we obtain 


_(^|0(w)) 

/C>£ 


^G|i;(w)) 


I-Pg) ^P£ 


{(t>n ( 

(l-G-f q)”'(g-q) 


3G) . (29) 


The above equation can be rewritten in a vector form 

P |^((a;)) = (T - Pg) PGQ |5G) . (30) 

Further on we use decoupling |0(a;)) = [P + Q]\G(}^)), 


Eq. (30), and Q\G{'^)) = a(a;)|S5G) to obtain 
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I^H) = (T-£ + q) ^ Q\g{u:)) =a{uj){l-C + Q) . 

We use Eq. (|26[) to determine the unknown function a{uj) 


(31) 


i^G) 


= — 


9G 

= —ILO ( 5G 


^C-C 


23S 


9g\ (3G|^(w)) +^(5G 


9/:-/: 


29E 


4>n ) {(j)n\Q{i^)) 


r 1 1 


^ ^ -1 ^ £| 



1 -f n - P£j P£ 


5G ) a(w). 


The final expression for the singular part of vector |0(a;)) then reads 

(3G) 


I^H) = 


-iujy/WP) l^G 


r ^ 11 



(i-£ + gj 


3G 


(T-£ + q) Vg) . 


(32) 


(33) 


Next, we rewrite the denominator in Eq. (33) in terms of the original two-particle vertex A. The expression in the 
square brackets can be written as 


S£kk' — £ 


kk' 


29Ek' 


— ■( Q^Akk' — - A 


ikk' 


3?Gk' 3?Ek' 


1 


+ 


2(3G2„)^„ 


+ 


9Gk(( 
5Gk5Gk^ 
{^Gi„ 


Ak'k"5iGk" — 2Ak'k"9Gk' 


— 3fiEk' ] + (^(^Akk"5iGk" — 2Akk''3'Gk' 


k"/k" 


(2 (sGk' 


A' 


k"k' 


-3Gk" 


k"k" 


- (9Gk"A, 


k"k' 


/SfJGk" 


k"k" 


S5Gk"3fiEk 


-5REk 3Gk 


5Gk' 


9Ek' ’ 


(34) 


where we employed a compact notation for momentum 
sums A^~^X)k/k = (/k)k- Expanding the inverse op¬ 
erator (1 — £ + Q) ^ in powers of £ — Q and using 
the properties Q|5G) =|5G), QC = Q, and the static 
Ward identity N~^ £kk'SEk' = SSk, we can explic- 

itly evaluate [ 1 — £ -I- Q] ^ I^^G) = lim„_>.oo £" |SG) = 
|3E) (5E|5G). We introduced a unit vector |3E) = 
(fV^(^)-i^^|k)3Ek in the last step. After a few 
additional manipulations we come to 


9G 


3 £-£ 


= ( 3G 


25E 


3 £-£ 


1 — £ -|- 
1 


-1 


3G 


25E 


3Ey (5E|3G) 


(SEkSGk)k(S^Gk)k 

2(5G2)^(5E2)^ ^ ^ 

and the final expression for the singular part of vertex 
t/(w) that has a rather simple form 





(36) 


two-particle behavior of the two-particle vertex becomes 
evident only in the g-dependent terms. 


V. DIFFUSION POLE AND STATIC 
DIFFUSION CONSTANT 

We use the singularity in the homogeneous vertex G{i-o) 
to determine the exact form of the diffusion pole and the 
diffusion constant controlling the low-energy asymptotics 
of the electron-hole correlation function. 


A. Derivation of the diffusion pole 

The uncorrelated part of the electron-hole correlation 
function <!>(£; w,q) is regular and hence the pole can 
emerge only in the vertex contribution 

^ E iE+)Gt (i^-)r^E 

kk' 

X G^,^{E+)G^,_{E.) (37) 


Notice that the low-frequency singularity of the homoge¬ 
neous two-particle vertex TkkK*^) is fully controlled by 
only one-particle functions. Equation ( |36|) is an alterna- 
tive version of the Velicky Ward identity^^ The genuine 


that will be found by summing the Bethe-Salpeter equa¬ 
tion over the fermionic momenta k and kk 

We first multiply Eq. (§ with 
AGk'(£; w,q)/AEk'(£;a;,q) and sum over k', which 
























































yields 


0«^(i?;a;,q) = (L. 


-k+k'' (^; q) 




AGk"(^;w,q) 


ASk"(i^;a;,q) 


+ G^,,iE+)GC,, {E.)g^,^{E;u;,ci) 


(38) 


J / k'' 


We introduced a reduced two-particle function 
g^^^iE-uj,^) = (5jf^k'(-^;‘^jq))k'- From now on, 

we will shorten AGk(i?;0,O) to AGk and similarly for 
ASk unless their full form is useful to indicate the lack 
of dependence on ui or q. 

In the next step we multiply Eq. ( [M| ) with 
AGk(-E; w, q), sum over k, and use the dynamical Ward 
identity AEk" (-E; w, q) = (E; w, q)AGk(E; w, q))j^ 

to eliminate the irreducible vertex We arrive at an 
integral equation for g^^{E\uj,q), 


AGk(E; q) - AEk(E; cc, q)G« (E+)G1 (E_) 


x0^^(£;;a;,q))^ = /AGk 


AEk(E;a;,q) 


(39) 


k \ ■■ ASk / k 

The product of two Green functions on the left-hand 
side is further rewritten in terms of differences, Eq. (12|, 
which yields 


/ [-UJ + Aek(q)] AGk(E^; ca, q) 
\ AEk(E; w, q) - w + Aek(q) 


g^^iE-co,q) 


= (^AGk 


AEk(E;a;,q) 

ASk 



(40) 


Since we are interested in the behavior of g^^{E;oj,q) 
only at small q and small w, it is sufficient to expand 
the factor standing at g^^{E-,(jj,q) in powers of q and 
w, and keep only terms up to q^ and w. Higher-order 
terms, including ujq, will be neglected. Because we expect 
g^^{E] w, q) to be singular at small q and uj, and because 
the expression on the right-hand side is regular, we can set 
Lo = 0 and q == 0 there right away. Furthermore, taking 
into account that Aek(q) is proportional to q, we can 
immediately write 


/ AGk(£^; 0,0) ^RA(p .,, N' 

■"^\ASk(E;0,O) 


Aek(q) 


AGk(E;;0,q) 


ASk(E;0,q) + Aek(q) 


g^^(E;a;,q) 


= (AGk)k- (41) 


We expand the multiplication factor at g^"^{E-,u},q) in 


Eq. (411 in powers of q = |q|. Since this function is already 


proportional to Aek(q)) we only need to expand Aek(q), 
AGk(E; 0, q) and ASk(E; 0, q) up to terms linear in q. 

The possibility to express everything in terms of real 
and imaginary parts of the retarded Green function and 
retarded self-energy comes from the fact that the argu- 
ment z in Eq. (13) is always z = E±i0. Inserting Eqs. (14) 


into the second term of Eq. (411 we arrive at 


Aek(q) AGk(E;0,q) 


= 9(q-Vk)|G^|2- 


AEk(£^;0,q) + Aek(q) 

V(q . Vk)|G«P [(q • Vk)3G« + 3(G^(q • Vk)E«)] , 


-ig^X>k 


(42) 


where we repeatedly used Eq. (12) at w = 0 and q = 0, 
that is, AGk/AEk = |GkP- Corning back to Eq. (41) we 
have 


-a;$(E;a;,q) +q((q- Vk)|G^P^^'^(i;;w,q))j^ 

-ig2(l?k^«^(E;cc,q))^ = (AGk)^, (43) 


where we introduced a reduced electron-hole correlation 
function 


^E-,u;,q) = {\G5\^g^^{E;u;,q))^ = 

(|G«pr«k^(£;;a;,q)|G^f)^^,. (44) 


The right-hand side of Eq. (43) can be evaluated with 
(AGk)^ = 2i($5Gk)j^ = —27rinF. For q = 0 we straight- 
forwardly find 


4>(a;) = (|G^P^^^(a;))k 


27rnF 

—iuj 


(45) 


For non-vanishing momentum q, we can rewrite Eq. (43) 


to 



iq{{q - Vk)\GQ\'^g^^^^,^{E;uJ,q})^^, + q^{Vkg+^^^,jE;u},q)) 


<+K_^ 


/kk' 




kk' 


^{E;uj,q) 


2'Kn-p 


(46) 
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where we used Eqs. (44) and (45) to represent <i)(i?;a;,q) 
in the denominator on the left-hand side. It does not 
depend on momentum q in the leading order. We now use 
the solution for vertices from the preceding section to 
obtain an explicit representation for the low-energy limit 
of the electron-hole correlation function. Using Eq. (20), 
the odd term in Eq. (46) reads 


9((q-Vk)|G^P^k+k; 

= g2^(q.Vk)|G«p 


(£;;w,q))kk, 

/k!fk”(q) 


l-C 

q ■ Vk» 

AEk» 


kk'k" 


q 

(47) 


which is also proportional to and contributes to the 
diffusion constant. 

We further extract the gradient of the one-electron 
functions from vertex £kk'(q)- We have 


(q) |G^ 


iV,£k+k;(g) = iV,Lk+K:^v 

^kk' 

-2!3G^W(G^,VS^,)] . 




Vk' 


(48) 


of consistent approximations for the diffusion constant 
needed to reach quantitative results. As a first step we 
choose the local mean-field approximation for which Ek = 
E and Akk' = SE/(3Gk")k" = A. Further on, £kk' = 
3Gk'/(5Gk")k"- Expression (51) for the static diffusion 
constant then reduces to the CPA resuit 


TrriFiA = ((q-Vk)^3G^)j^ . (52) 

This formula holds for any approximation with a local 
vertex A. 

Approximations with the local irreducible vertex do 
not improve upon the Drude or CPA diffusion constant. 
Only a non-local irreducible vertex can make a difference. 
Since it is complicated to evaluate the full momentum 
dependence of the self-energy resulting from momentum 
dependence of the vertex function, we consider only the 
momentum dependent self-energy as a correction to the 
local term. We denote 

E^iE) = E«(E) + SE^iE) , (53a) 

A«^(E;a;,q) = X^^iE;u;) + SA5^,{E-,LU,q) . (53b) 

and take the non-local correction from the maximally 
crossed diagrams, that is 


Since vertex £kk'(q) depends only on even powers of 
momentum g, V^Lkk' (q) = 0. With this representation 
we obtain the singular part of the low-energy limit of the 
electron-hole correlation function in form of the canonical 
diffusion pole 


$(E;a;,q) 


27rnF 

—iw -I- D{uj)q^ 


(49) 


with the diffusion constant 


D{u:) = 


lUJ 

27Tni 


(q - vk)|G| 


R\2 


l-C 


[9G^,q-Vk' 

kk' 


+5 (G^,q-Vk'E^,)]tI«^(..)) 


kk'k'' 


(50) 


Equations ( |4^ and ( [50| determine the exact canonical 
form of the diffusion pole in the electron-hole (density- 
density) correlation function. 


B. Static diffusion constant 


dA«fk'_, [E; cu, q) ee SA{E; uj,k + k') 

^ X^^{E;u;rx^{E;oj,k + k') 

l-XRA{E;oj)x+{E;uj,k + k') ■ ^ > 

We used the reduced electron-electron bubble 


X^^iE-,to,q)=x^^{E;uj,q) 

-G^(E+)G^(E_) = ^5]G«^(E+)GV('E;-) 

k 

-G^(E+)G^(E_) . (55) 


The non-local correction to the CPA irreducible vertex 
leads to a correction to the CPA the self-energy. 


3dEk(E) = ^ 

3fidEk(E) = P 


^MS()(E)3Gk'(E) , 



duj S5(5Ek(a;) 
TT U) — E 


(56a) 

(56b) 


An explicit representation for the static diffusion con¬ 
stant D — D{0), that is, the explicit limit o; —>■ 0 in 
Eq. (50) is obtained from Eq. (50) when we use the sin¬ 
gular part of vertex ^^"^(w) from Eq. (36). The static 
diffusion constant is then expressed via a Kubo-like for¬ 
mula with the full two-particle vertex 


nupE = ((q • Vk)|GkP [lV(5k.k' + rkk'|Gk'P] 

X [9Gk'q • Vk' + 3 (Gk'q • Vk'Ek')] 3Sk')kk' > (51) 

where we used [1 — £]kk, = A^<5k,k' + rkk'|Gk'p. This 
exact expression is the starting point for the derivation 


Generally, this correction to the mean field is justified for 
weak disorder where the momentum dependence of the 
self-energy is expected to be small. The one-electron prop- 
agators contain only the local, unperturbed self-energy. 

The non-local irreducible vertex (5Ak^ depends only 
on the sum k -f k' and contains the Cooper pole at k -f 
k' = 0. The contribution from this pole to the diffusion 
constant will be dominant in low dimensions. We use this 
dominance and also neglect the momentum dependence 
of the self-energy (the correction i5Ek). To determine 
the diffusion constant, we have to evaluate the following 
expression 
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Wk 


1 - A |G|"j = vInS^m' + E E E • • • E ^kA(k + ki)|Gk, |"A(ki + k2)|Gk, |^ .. <5k„,k' 




Z)k' 


k, ko 


oo 

= w^A^4.k' + E E E ■ • • E^kA(qi)|Gq,_kpA(q 2 )|Gq,_qi+kP ... A(q„)|Gk„p4„,k'Wk' (57) 


n—1 


qi q 2 


where k„ = + (—l)"k. Since the vertex 

correction (5A(q) contains the Cooper pole, the doniinant 
contribution to the integrals in Eq. ( |57[ ) comes from the 
values of non-singular functions at qi = 0. The multiple 
sums (integrals in the thermodynamic limit) then factorize 
and we obtain a renormalized diffusion constant in a closed 
form 


TTUpD = 


(q-Vk)^9G2 


l + |GkP(M«^(g)) 


(58) 


'i/ k 


where 

(^A«^('7))q 


xr(o) 


(0) - |GP 

X?^(0)-X?^(q) 


- 1 


■ (59) 


The mean-field self-energy E^(G) is determined from the 
Soven equation 


1 = 


1 


1 + [E(E) - E,] G{E ), 
and the corresponding irreducible vertex is 
3E(A;) 1 1 


(60) 




^G{E) (|Gk(E)P)k |G(E)P 


[l -/ 

1 

yi 

[ \ 

1 + (S(F;) - V,) G{E) 

/ av 


(61) 


The expression for the static diffusion constant in 
Eq. (58) is similar to the expression for vertex correc- 
tions to the mean-field conductivity due to maximally 
crossed diagrams derived in Ref. [231 Although they both 
lead to a non-negative resuit and correctly describe weak 
localization in one- and two-dimensional Systems, they are 
strictly justified only in the weak-scattering limit where 
the non-local corrections are negligible compared to the 
CPA self-energy. The diffusion constant from Eq. (58) 
is nevertheless attractive in that it offers an easily ac- 
cessible qualitative assessment of the impact of vertex 
corrections to the mean-field, Drude resuit without the 
necessity to go beyond the Soven equation for the local 
self-energy. Moreover, this resuit can easily be extended 
to multi-orbital models and used in realistic calculations. 


As an example, we calculated the static diffusion con¬ 
stant for a three-dimensional cubic lattice with the ran- 
dom potential corresponding to a binary alloy, and with 
a potential given by a box distribution. We compared the 
CPA resuit, Eq. (52), with a first-order vertex correction 
from the denominator on right-hand side of in Eq. (58), 
and the full solution of Eq. (58). The results are plotted 
in Eig. We can see that the values of the diffusion 
constant do not differ much deep inside the conduction 
bands, but the vertex corrections make the Drude resuit 
unstable near the band edges. Full solution of Eq. (58) 
regularizes the correction so that the diffusion constant re- 
mains nonnegative. Neither of these Solutions is, however, 
capable to describe the mobility edge or the Anderson 
localization transition. The diffusion constant decreases 
with increasing disorder but it never reaches zero for fi¬ 
nite values of the disorder strength in this approximation, 
see Fig. One needs a self-consistent approximation for 
the two-particle vertex in order to describe vanishing of 
diffusion at the Anderson localization transition in three 
and higher spatial dimensions. 


VI. CONCLUSIONS 

Charge diffusion is a macroscopic phenomenon charac- 
terized by a specific form of the low-energy limit of the 
density-density response function. This function displays 
a diffusion pole where the long-range spatial fluctuations 
are controlled by a diffusion constant. Unfortunately, the 
density response function and the diffusion constant are 
not elementary objects of the microscopic theory of quan¬ 
tum diffusion. They are composite objects represented by 
the sum over the fermionic momenta of the two-particle 
Green function. This fact makes the full description of 
quantum diffusion from first, quantum-mechanical princi¬ 
pies a demanding task that is very difhcult to accomplish 
without inconsistencies and ad hoc or unjustified steps. 
The objective of this paper was to derive an exact repre- 
sentation of the diffusion pole and the diffusion constant 
in a consistent way from the elementary objects of the 
diagrammatic perturbation theory. 

A ubiquitous severe problem of the perturbation the¬ 
ory of non-interacting quantum random systems is the 
inability to obey the Ward identity when we go beyond 
the local, mean-field approximation. Whatever approach 
we choose, either an approximation for the self-energy 
(one-particle approach), or an approximation for the irre- 
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FIG. 1: (Color Online) Diffusion constant for a binary alloy with component concentrations 30% and 70%, and with the difference 
between the local potentials AV = Wo, where 2Wo is the band with of the lattice without disorder (left panel). Right panel 
shows analogous calculation for a random potential corresponding to the box distribution of width 2Wo. Plotted are the values 
of the diffusion constant for the CPA, Eq. (52 I (dotted, red), CPA with the first-order vertex correction (dashed, green), and the 
full solution of Eq. (|58|l (solid, bine) on a simple cubic lattice. 



E/Wq 

FIG. 2: (Color Online) Diffusion constant from Eq. ( |58| l for 
various box distributions of width 2Wo,4:Wo, and 81Fo on 
a simple cubic lattice. The diffusion constant decreases with 
increasing disorder, but it does not reach zero for finite disorder 
strengths. 


ducible vertex (two-particle approach), the Ward ideiitity 
cannot be guaranteed and is only used as a consistency 
check. In this paper, we made a substantial step forward 
and clarified the way the two-particle irreducible vertices 
from the perturbation theory should be treated in or- 
der to restore the full dynamical Ward identity and all 
macroscopic conservation laws of measurable quantities. 
We used the Bethe-Salpeter equation with the irreducible 
vertex complying with the Ward identity and derived 
an exact low-energy singular asymptotics of the full two- 
particle vertex. We then represented the diffusion pole 


and the diffusion constant in the low-energy asymptotics 
of the electron-hole correlation function via the func- 
tions resulting from the perturbation expansion, that is, 
the self-energy and the two-particle irreducible vertex. 
It follows from our analysis that the singular structure 
of the two-particle vertex is more complex than that of 
the electron-hole correlation function. Consequently, the 
diffusion constant cannot be directly pulled into the per¬ 
turbation theory as a parameter of a two-particle function 
to be self-consistently determined. 

The diffusion constant remains a parameter set apart 
of the perturbation theory. It contains only a reduced 
Information from the two-particle vertex. It is represented 
as a matrix element of the full two-particle vertex via 
a Kubo-like formula. It follows from the exact represen- 
tation derived here that any local approximation to the 
two-particle vertex irreducible in the electron-hole scat- 
tering channel results in the mean-field diffusion free of 
vertex corrections. Only non-local approximations with 
momentum-dependent irreducible vertices produce cor¬ 
rections to the Drude term. The contribution from the 
maximally crossed diagrams was explicitly calculated as- 
suming dominance of the low-energy singularity in the 
two-particle vertex. We obtained a closed form of the 
static diffusion constant containing the leading vertex 
corrections and describing the weak localization that can 
easily be generalized to multi-orbital models and can be 
used in realistic calculations. Since the diffusion constant 
is not an integral part of the perturbation theory, finding 
a criterion of Anderson localization in terms of the irre¬ 
ducible vertex functions is stili an open problem. It also 
remains unclear what is the minimal approximation be- 
ing able to describe qualitatively correctly the Anderson 
localization transition. A framework for addressing these 
























12 


questions was set in this paper. 
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